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ABSTRACT 

The  popularity  of  micropropulsion  system  development  has  led  to  renewed  interest  in  the  determination  of 
propulsive  properties  of  orifice  flows  since  micronozzle  expansions  may  suffer  high  viscous  losses  at  low 
pressure  operation.  The  mass  flow  and  relative  thrust  for  an  under  expanded  orifice  is  measured  as  a 
function  of  orifice  stagnation  pressure  from  0.1  to  3.5  Torr.  Nitrogen,  argon,  and  helium  propellant  gases 
are  passed  through  a  1.0  mm  diameter  orifice  with  a  wall  thickness  of  0.015  mm  .  Near- free  molecule, 
transitional  and  continuum  flow  regimes  are  studied.  The  relative  thrust  is  determined  by  a  novel  thrust 
stand  designed  primarily  for  low  operating  pressure,  micropropulsion  systems.  It  is  shown  that  the  thrust 
indications  obtained  from  the  stand  are  a  function  of  the  facility  background  pressure,  and  corrections  are 
made  to  determine  the  indicated  thrust  for  a  zero  background  pressure  with  nitrogen  as  propellant.  Highly 
repeatable  (within  1  %)  indicated  thrust  measurements  are  obtained  in  the  thrust  range  from  5  to  500  pN. 


Nomenclature 

a  -  speed  of  sound  (m/s ec) 

A  -  area  (m2) 

c(y)  -  constant  dependent  on  ratio  of  specific 
heats 

c’  -  mean  molecular  thermal  speed  (m/sec) 
Cd  -  discharge  coefficient 
d  -  diameter  (m)  ' 

g  -  gravitational  constant  (=  9.8  m/sec2) 

Isp  -  specific  impulse  (sec) 
k  -  Boltzmann’s  constant  (=  1.38  x  10‘23  J/K) 
Kn  -  Knudsen  Number 
m  -  mass  (kg) 

-  mass  flow  (kg/sec) 
n  -  number  density  (m"3) 
p  -  pressure  (Pa) 
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Re  -  Reynolds  Number 
rp  -  radius  of  penetration  (m) 
t  -  thickness  (m) 

T  -  temperature  (K) 

a  -  orifice  transmission  probability 
y  -  ratio  of  specific  heats 
X  -  mean  free  path  (m) 
p  -  viscosity  (Ns/m2) 
p  -  mass  density  (kg/m3) 

A  -  thrust  stand  deflection  (arb.  units) 

subscripts 

b  -  facility  background 
fm  -  free  molecule 
L  -  limit  (theoretical  maximum) 
meas  -  experimentally  measured 


o  -  stagnation  region 
p  -  plenum 

t  -  orifice  or  nozzle  throat  property 
theor  -  theoretical 

superscripts 

*  -  orifice  plane  (sonic)  region 


1.  Introduction 

In  recent  years,  micropropulsion  systems  have 
been  developed  to  address  the  need  for  highly 
mobile  micro-  and  nanospacecraft  A  wide  array 
of  concepts  will  require  the  expansion  of 
propellant  gases  through  microscale  geometries 
(e.g.  micronozzles).  Because  of  the  volume  and 
power  restrictions  associated  with  storing  or 
producing  high  pressures  on  micospacecraft, 
many  micropropulsion  systems  will  operate  at 
relatively  low  pressures  in  the  transitional  flow 
regime.1  The  Reynolds  number  gives  a  measure 
of  the  flow  efficiency  in  terms  of  viscous  losses. 
The  Reynolds  number  at  a  nozzle  throat  or  an 
orifice  is  given  by 


* 

Re 


*  *  j 

p  a  dt 
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(1) 


Lower  Reynolds  number  implies  higher  viscous 
flow  losses.  Microspacecraft  propulsion  systems 
may  inherently  operate  in  low  Reynolds  number 
regions  due  to  relatively  low  operating  pressures 
and  small  characteristic  dimensions. 


Figure  1  shows  the  specific  impulse  as  a  function 
of  distance  through  a  conical  micronozzle 
geometry  with  a  throat  diameter  of  dt  =  27.7  /mi. 
Navier-Stokes  and  Direct  Simulation  Monte 
Carlo  numerical  simulations  have  been 
performed  at  two  different  stagnation  pressures 
p0  =  106  Pa  and  105  Pa.2  For  cold  gas  operation 
(T0  =  300  K),  the  corresponding  Reynolds 
numbers  are  1300  and  130  respectively.  As  Fig. 
1  shows,  the  specific  impulse  at  the  nozzle  exit  is 
approximately  14%  higher  than  at  the  nozzle 
throat  for  Re*  =  1300.  However,  the  specific 
impulse  at  the  exit  is  only  about  5%  higher  than 
at  the  nozzle  throat  for  Re*  =130.  The  reduction 
of  efficiency  in  the  micronozzle  geometry  as  the 
Reynolds  number  decreases  is  due  to  viscous 
losses  near  the  nozzle  walls. 


There  is  currently  a  large  effort  being  devoted  to 
the  fabrication  of  micronozzles  with  throat 
diameters  on  the  order  of  one  to  tens  of 
micrometers  (microns).3-5  As  Fig.  1  indicates, 
weighing  the  increase  in  performance  versus 
fabrication  complexity.  At  low  Reynolds 
numbers,  the  micronozzle  geometry  may  in  fact 
degrade  performance,  and  expansion  of 
propellant  from  a  simple  thin  walled  orifice  may 
be  a  good  compromise  between  efficiency  and 
system  complexity.  Some  microelectro¬ 
mechanical  systems  (MHMS)  fabricated  nozzle 
geometries  involve  planar  or  rectangular  throats 
(i.e.  not  axisymmetric).3  Numerical  studies  have 
shown  that  flows  generated  near  the  side  walls 
can  result  in  even  higher  inefficiencies.6 

The  flow  complexities  from  sonic  orifices  have 
been  studied  for  several  years.7-11  However,  the 
determination  of  the  thrust  generated  from  gas 
expanding  through  an  orifice  is  an  area  that  has 
received  little  attention.  The  advent  of 
micropropulsion  systems  has  renewed  interest  in 
the  determination  of  propulsive  properties  of 
orifice  flows  since  micronozzle  expansions 
appear  to  have  major  viscous  losses. 

This  manuscript  explores  the  thrust  generated  by 
an  orifice  expansion  at  relatively  low  Reynolds 
number.  Besides  a  propulsion  system  in  its  own 
right,  these  orifices  are  also  being  investigated  as 
a  reliable  means  of  calibrating  micro-Newton 
thrust  stands.12 


2.  Theory 


To  assess  the  performance  of  the  orifice 
expansion  in  terms  of  propulsive  parameters, 
properties  at  the  entrance  plane  of  the  orifice  are 
calculated  from  known  stagnation  values.  The 
ratios  of  pressure,  density,  temperature  and 
velocity  for  inviscid  flow  are 
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The  theoretical,  inviscid  flow  value  for  the 
orifice  mass  flow  is 

A?=  p*  a  At  (6) 


Viscous  effects  can  be  measured  in  terms  of  a 
discharge  coefficient  defined  by 
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where  Mlheor  is  calculated  from  Eq.  (6).  The 

theoretical  thrust  produced  by  the  orifice  is  then 
given  by 


3.=  itithenrO*  +P*At  (8) 
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The  constant  c(y)  is  equal  to  1.16  and  1.14  for  y 
=  1.4  and  1.67  respectively.  Equation  (8)  is  only 
valid  for  the  relatively  low  stagnation  pressure 
range  under  investigation  in  this  study  where  the 
effects  of  the  external  flow  expansion  can  be 
neglected.  Therefore,  it  is  expected  that  the 
thrust  produced  by  the  orifice  gas  flow  is 
relatively  independent  of  the  propellant. 

A  measure  of  the  propulsive  efficiency  is  given 
by  the  specific  impulse  as 
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For  free  molecule  flow,  the  Knudsen  number 
defined  by 


Knt  = 


4a 

dt 
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is  relatively  high  (Kn  >  10).  This  is  realized  at 
very  low  stagnation  pressure  operation  where  the 
molecule  mean  free  path  is  larger  than  the  orifice 


diameter.  The  free  molecule  mass  flow,  thrust 
and  specific  impulse  are  given  by 
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For  the  thin  walled  orifice  used  in  this  study,  the 
theoretical  transmission  probability  through  the 
orifice  a  is  near  unity.11 

3.  Experiment 


The  orifice  used  in  this  study  has  a  diameter  of  dt 
=  1.0  mm  and  a  thickness  of  tt  =  0.015  mm 
giving  a  t/d  =  0.015.  For  t/d  =  0.015,  the 
transmission  probability  (a  in  Eqs.  (12)  and  (13)) 
is  very  close  to  unity.13  The  orifice  is  machined 
by  conventional  means  in  a  Tantalum  shim 
which  is  attached  to  an  aluminum  plenum  as 
shown  in  Fig.  2.  The  aluminum  plenums  are 
attached  to  a  torsional  thrust  stand  shown  in  Fig 
3.  The  thrust  measurements  involve  sensing  the 
angular  displacement  resulting  from  a  torque 
(thrust  force  x  radial  distance)  applied  to  a 
damped  rotary  system.  The  present  method  for 
detecting  angular  deflection  is  to  measure  the 
linear  displacement  at  a  known  radial  distance 
using  a  Macro  Sensors,  Inc.  linear  voltage 
differential  transducer  (LVDT).  The  total  linear 
movement  of  the  arm  is  approximately  0.5  mm 
for  a  2mN  thrust  level  which  corresponds  to  less 
than  0.1°  angular  deflection.  The  detailed 
operational  characteristics  of  this  thrust  stand  is 
the  topic  of  earlier  work.7 

The  thrust  stand  is  placed  inside  the  CHAFF-II 
facility,  a  steel  vacuum  chamber  pumped  by  a 
Roots  blower  system  with  a  pumping  speed  of 
2000  L/sec  for  nitrogen.  Ultimate  pressures 
achievable  in  CHAFF JI  are  approximately  1.0  x 
10*4  Torr. 

The  propellant  is  introduced  into  the  orifice 
plenum  through  an  adjustable  needle  valve 
located  downstream  of  an  MKS®  mass  flow 
meter.  In  the  experimental  configuration,  the 
mass  flow  meter  operated  in  the  continuum 
regime  through  the  pressure  range  studied. 
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Nitrogen,  argon  and  helium  are  used  as 
propellant  gases  in  this  study. 

4.  Results 

Figure  4  shows  the  discharge  coefficient  as  a 
function  of  the  Reynolds  number  for  nitrogen 
and  argon,  At  lower  Reynolds  number,  the 
measured  values  should  asymptotically  approach 
the  theoretical  free  molecule  limit  of  0.583  and 
0.549  for  y  =  1.4  and  1.67  respectively.  For  the 
range  of  Reynolds  number  shown  in  Fig.  4,  the 
orifice  stagnation  pressure  ranges  from  0.1  to  3.5 
Torr.  The  discharge  coefficient  asymptotes 
towards  0.9  for  higher  Reynolds  number  flows. 
The  maximum  discharge  coefficient  for  nitrogen 
is  consistent  with  previous  studies  which  show 
that  for  sharp  edged  orifices  (t/d  -  0.01),  the 
maximum  discharge  coefficient  is  approximately 
0.9  for  Re*  upwards  of  104.10 

Figure  5  shows  typical  thrust  stand  traces  from 
the  LVDT.  The  traces  are  for  nitrogen  at  p0  = 
0.1  and  3.5  Torr  respectively.  Using  Eq.  (8)  with 
the  measured  mass  flow  and  stagnation  pressure, 
the  theoretical  range  of  thrust  shown  in  Fig.  5  is 
from  approximately  7.9  fiN  (p0  =  0.1  Torr)  to 
430  jttN  (p0  =  3.5  Torr).  As  indicated  by  the 
traces,  the  signal  voltage  at  p0  =  0.1  Torr  is 
approximately  15  to  20  times  above  the  noise.  In 
this  study,  the  noise  is  generally  dominated  by 
the  minimum  resolution  of  the  digitizer  (i.e.  bit 
noise)  and  not  an  intrinsic  property  of  the  thrust 
stand.  This  suggests  that  the  thrust  stand  should 
be  capable  of  measuring  thrust  levels  below 
1.0  fxN  with  minor  modifications. 

The  indicated  thrust  measurements  obtained 
from  the  thrust  stand  have  a  standard  deviation 
of  6.65  x  10'3  and  1.04  x  10*2  for  five  runs  on 
nitrogen  and  argon  at  p0  =  1.0  Torr  respectively. 
Based  on  the  mean  deflection  of  the  stand,  this 
translates  to  a  statistical  error  of  0.59%  for 
nitrogen  and  0.95%  for  argon.  In  all,  the 
standard  deviation  is  within  1%  of  the  mean 
deflections  over  the  range  of  stagnation  pressures 
from  0.1  to  3.5  Torr. 

The  measured  linear  deflection  from  the  thrust 
stand  mounted  LVDT  is  shown  as  a  function  of 
orifice  stagnation  pressure  for  nitrogen 
propellant  in  Fig.  6.  Similar  plots  are  shown  for 
argon  and  helium  in  Figs.  7  and  8  respectively. 
Figure  9  shows  the  nitrogen  deflection  as  a 
function  of  mass  flow. 


Figure  10  shows  the  thrust  stand  deflection  for 
both  nitrogen  and  argon.  Because  the  range  of 
operating  Reynolds  number  for  argon  and 
nitrogen  are  similar,  the  results  overlap  as 
expected  from  Eq.  (9).  However,  the  level  of 
agreement  only  holds  for  propellants  which  have 
a  similar  Reynolds  number  for  a  given  stagnation 
pressure  (i.e.  similar  viscous  effects  as  a  function 
of  p0). 

Because  the  orifice  Reynolds  number  (Eq.  (1))  is 
a  factor  of  three  lower  for  helium  than  nitrogen, 
it  is  expected  that  the  helium  thrust  (deflection) 
level  would  be  somewhat  lower  at  a  given 
stagnation  pressure.  The  helium  results  are 
compared  to  nitrogen  deflections  in  Fig.  11.  As 
expected,  the  deflection  is  lower  than  that  for 
nitrogen  propellant  at  similar  stagnation 
pressures  due  to  higher  viscous  losses. 

It  is  known  that  the  deflection  for  a  given  orifice 
stagnation  pressure  is  dependent  on  the 
background  pressure  of  the  facility.  Figure  12 
shows  the  measured  deflection  for  given  orifice 
stagnation  pressures  as  a  function  of  the  chamber 
background  pressure.  The  absolute  deflection 
for  a  given  stagnation  pressure  decreases  as  the 
background  pressure  increases.  The  slope  and 
intercept  of  the  data  is  used  to  correct  the  results 
in  Fig.  6  to  a  zero  background  pressure  in  the 
following  section. 

5.  Discussion 

Indications  of  Rarefied  Flow 

Figure  13  shows  the  effects  of  different  flow 
regimes  on  the  measured  thrust  (deflection)  from 
a  sharp-edged  orifice.  At  high  Knudsen  number, 
the  measured  thrust  follows  the  free  molecule 
theory  solution  (Eq.  (13)).  As  the  transitional 
flow  regime  is  reached,  the  measured  flux 
deviates  from  free  molecule  theory  and  begins  to 
asymptote  to  the  inviscid  flow  theory  solution 
(Eq.  (8))  as  the  continuum  regime  is  reached. 
Because  the  discharge  coefficient  only  reaches  a 
maximum  of  approximately  0.9  for  sharp-edged 
orifices  (Fig.  4  and  Ref.  10),  the  inviscid  flow 
solution  is  never  quite  met. 

For  helium,  the  flow  Knudsen  number  at  p0  =  0.1 
Torr  is  approximately  1.4.  The  helium 
deflection  data  in  Fig.  14  shows  the  theoretical 
lines  for  free  molecule  and  inviscid  continuum 
thrust.  The  measured  data  is  bounded  by  the  free 
molecule  and  continuum  solutions.  At  the  lower 
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operating  pressure,  the  data  follows  the  free 
molecular  slope  and  tends  toward  the  inviscid 
continuum  solution  at  higher  pressures.  The  data 
does  not  reach  the  inviscid  solution  at  the 
maximum  Reynolds  number  (p0  =  3.5  Torr)  of 
about  27. 

Figure  15  shows  a  similar  result  for  nitrogen 
propellant  flows.  The  maximum  Knudsen 
number  reached  is  approximately  0.5  at  p0  =  0.1 
Torr.  Therefore,  only  a  slight  indication  of  the 
data  trending  towards  the  rarefied  solution  is 
expected  as  shown  in  Fig.  15  for  p0  <  0.5  Torr. 
At  the  higher  operating  pressures,  the  nitrogen 
indicated  thrust  closely  follows  the  inviscid 
continuum  solution.  This  is  consistent  with  the 
discharge  coefficient  results  shown  in  Fig.  4 
which  indicates  a  discharge  coefficient  near  0.85 
for  Re*  >  40. 

Facility  Background  Pressure  Corrections 

Because  the  facility  background  pressure  is  made 
up  of  two  components  (laboratory  air  and 
propellant),  corrections  to  the  thrust  stand 
deflections  can  only  be  approximated  for  the 
nitrogen  flow  cases  where  the  propellant  and  the 
laboratory  backgrounds  are  similar. 

The  mechanism  for  lower  thrust  deflection  as  a 
function  of  increased  background  pressure  is 
shown  schematically  in  Fig.  16.  For  no  orifice 
flow,  the  background  pressure  exerts  an  equal 
force  on  the  front  and  back  sides  of  the  orifice 
plenum  (equilibrium).  As  flow  is  introduced 
through  the  orifice,  the  resulting  jets  acts  like  an 
ejector  pump  similar  to  that  of  the  oil  vapor  in  a 
vacuum  diffusion  pump.13,14  Collisional 
“removal”  of  the  background  gas  by  the  orifice 
plume  results  in  a  lower  background  pressure  on 
the  jet  side  of  the  plenum  then  on  the  back  side. 
The  pressure  difference  exerts  a  force  on  the 
orifice  in  a  direction  opposite  of  the  thrust  vector 
produced  by  the  jet.  Since  the  gas  density  in  the 
plume  is  relatively  high  compared  to  the 
background  gas  density  in  the  vicinity  of  the 
orifice,  the  source  flow  can  effectively  prevent 
background  molecules  from  penetrating  the 
orifice  plume  and  striking  the  front  surface  of  the 
orifice  plenum  for  reasonable  background 
pressures  (pb  <  1  x  10'3  Torr).1:>  This  suggests 
that  the  deflection  as  a  function  of  the 
background  pressure  should  be  linear  as  a  first 
order  approximation  as  Fig,  12  indicates,  or 


dA 

- =  constant 

dPb 

(Vi,=o)-  (V6  )=  PbAeJf  (15) 

The  effective  area,  Aeff,  is  a  measure  of  the  area 
being  utilized  to  cause  the  deflection  opposite  to 
the  thrust  vector.  Assuming  that  the  background 
pressure  in  front  of  the  orifice  plenum  is  zero, 
the  force  exerted  on  the  back  side  of  the  orifice 
plenum  (Ap  =  22.75  cm2)  is  approximately  100 
/xN  for  pfa  =  3.3  x  10‘4  Torr.  Therefore,  the 
background  pressure  “negative  thrust”  effect  can 
be  a  significant  fraction  of  total  “positive  thrust” 
produced  by  the  orifice. 

Figure  17  shows  the  nitrogen  thrust  stand 
deflection  for  a  corrected  zero  background 
pressure.  This  data  is  derived  from  the  slopes  of 
the  deflection  versus  background  pressure  curves 
in  Fig.  12  and  the  data  obtained  in  Fig.  6. 

6.  Conclusions 

The  maximum  of  the  discharge  coefficient  in  this 
study  was  approximately  0.86  indicating  that 
viscous  effects  are  present  throughout  the 
operational  pressure  range.  The  momentum  flux 
measurements  from  a  sharp-edged  orifice  (t/d  = 
0.015)  show  a  linear  dependence  on  the 
stagnation  pressure  and  mass  flow  for  Re*  >  20. 
In  this  Reynolds  number  range,  the  data 
approaches  the  inviscid  flow  solution.  Below 
Re*  =  20,  rarefied  effects  begin  to  influence  the 
thrust  results,  and  the  data  approaches  the  free 
molecule  theoretical  solution  as  seen  in  the 
deflection  data  for  helium  at  low  stagnation 
pressure. 

The  thrust  stand  deflection  is  a  function  of  the 
facility  background  pressure.  Corrections  to  the 
deflection  data  must  be  made  by  extrapolating  to 
a  zero  background  pressure  condition.  The 
correction  is  upwards  of  about  40%.  Operating 
the  thrust  stand  with  at  least  an  order  of 
magnitude  lower  ultimate  facility  pressure  is 
desired. 

Highly  repeatable  (within  1  %)  indicated  thrust 
measurements  are  obtained  in  the  thrust  range 
from  5  to  500  /xN.  There  is  also  data  to  suggest 
that  the  thrust  stand  will  be  capable  of  measuring 
thrust  levels  below  1.0  /xN  with  minor 
modifications. 
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An  orifice  which  achieves  discharge  coefficients 
near  unity  can  be  used  as  a  calibration  source  for 
thrust  stand  measurements.  An  orifice  has  the 
advantage  of  being  able  to  obtain  in-situ 
calibration  without  moving  parts.  The  accuracy 
of  the  results  may  be  much  higher  than 
traditional  means  of  thrust  stand  calibration  since 
the  inviscid  flow  solutions  are  expected  to  be 
accurate  for  CD  ~  1  and  the  flow  can  be 
simulated  numerically  as  well. 
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Fig.  1:  Specific  impulse  along  the  axis  of  a  MEMS  fabricated  micronozzle,  (ref.  1) 


Fig.  2:  Orifice  geometry  (t/d  =  0.015). 
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Fig.  3:  Micro-Newton  thrust  stand  installed  in  CHAFF- II  vacuum  facility. 
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Fig.  4:  Discharge  coefficient  as  a  function  of  orifice  Reynolds  number  for  argon  and  nitrogen  gas  flows 
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Fig.  5:  LVDT  signal  from  thrust  stand.  (A)  Nitrogen  flow,  p0  =  0.1  Torr;  (B)  Nitrogen  flow,  p0  =  3 
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Fig.  6:  Deflection  as  a  function  of  stagnation  pressure  for  nitrogen. 
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Fig.  7:  Deflection  as  a  function  of  stagnation  pressure  for  argon. 
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Fig.  10:  Deflection  as  a  function  of  stagnation  pressure  for  nitrogen  and  argon  gas  flows. 
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Fio.  pi;  Deflection  as  a  function  of  stagnation  pressure  for  nitrogen  and  helium. 
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Fig.  12:  Deflection  as  a  function  of  chamber  background  pressure  for  nitrogen  gas  flows. 
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Fig.  13:  Illustration  of  flow  regime  effects  on  the  measured  relative  thrust. 
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Fig.  14:  Deflection  as  a  function  of  stagnation  pressure  for  helium. 


0.0  0.50  1.0  1.5  2.0  2.5  3.0  3.5  4.0 

Stagnation  Pressure  (Torr) 


Fig.  15:  Deflection  as  a  function  of  stagnation  pressure  for  nitrogen  gas  flow. 
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Fig.  17:  Corrected  deflection  to  zero  background  pressure  for  nitrogen. 
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